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ABSTRACT 


In a computer simulation experiment, acoustic tomography is assessed as a 
means of measuring the seasonal flux of heat advected by the Norwegian 
Atlantic Current. Oceanic heat flux has traditionally been measured by various 
direct or indirect techniques that are prone to error or large uncertainty. The 
tomographic technique offers distinct advantages over conventional methods in 
that temperature and current fields, that combine to yield heat flux in the ocean, 
can be determined at various spatial and temporal scales. The adequacy of the 
tomographic technique thus hinges on the question of "how well can the 
temperature and current be resolved spatially?" The spatial resolution of 
tomography varies with array size, number of transceivers used and the 


characteristics of the sound channel. 


In the assessment we use the General Digital Environmental Model 
(GDEM), a climatological data base, to simulate an ocean area 550 x 550 km2 off 
the Norwegian Coast. Resolution and variance analysis are performed on two 
circular arrays consisting of six transceivers. An important finding is that the 
horizontal resolution lengths of the current and temperature fields differ. For a six 
element array the horizontal resolution length is approximately one fifth the array 
diameter for the current field , whereas for the temperature field it is one sixth the 


array diameter. 


We then generate synthetic travel time data that have embedded within them 
temperature and current signals as well as random noise. We invert the synthetic 
travel time data to form estimates of the original fields using a linear optimal 


estimator based on the Gauss-Markoff theorem. We relate the sound speed 


perturbation field to potential temperature and compare these estimates to the 
original values. Finally we use the estimated fields to compute heat flux across a 
transect located within the array. We compare the actual to the estimated heat 
flux to asses the quality of the tomographically derived value. We have found that 
the quality of the heat flux estimates depends critically on how well the flow field 
is resolved. A six element array can adequately resolve the current in the 
Norwegian Sea, provided that its diameter is shorter than 250 km. Such an array 
is able to measure net heat flux through a transect at the center of the array with 


only a 10% error. 
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I. INTRODUCTION 


Ocean acoustic tomography is a technique that uses a fixed array of 
transceivers to measure travel times of pulses along acoustic multi-paths. These 
measurements contain the information necessary to allow determination of 
temperature t(x,y,z,t) and current V(x,y,z,t) structure. The concept was first 
introduced by Munk and Wunsch (1979) as a means of large scale monitoring of 
mesoscale processes occurring in the ocean. Since that time acoustic 
tomography has become a viable tool for measuring the oceans physical 
properties and has greatly enhanced our ability to probe beneath the ocean 
surface. Ocean acoustic tomography has been successfully demonstrated in a 
variety of applications such as the mapping of a mesoscale eddy field by 
Cornuelle et al (1983), planetary wave analysis by Chiu and Desaubies (1987) 
and surface wave spectra analysis by Lynch et al (1987). This thesis explores 
yet another application of ocean acoustic tomography, namely the measurement 


of ocean advective heat flux. 


A. PURPOSE 

The world's oceans play a major role in the redistribution of energy to 
maintain the earth in thermal balance. The equatorial regions of the earth receive 
a net excess of energy from insolation while the poles lose more than they 
receive. Excess energy at the equator that cannot be reradiated back into space 
must be advected away by the atmosphere and ocean to maintain a local thermal 
equilibrium. The net energy flux is poleward over the entire globe. Many attempts 


have been made over the past four decades to accurately compute ocean heat 


flux. According to Vonder Haar and Oort (1973), the oceans are responsible for 
up to 40% of the heat flux. However past measurements of the oceans 
contribution have been in error by as much as 70%. The poleward flux of heat has 
a governing effect on several processes including, ocean circulation, weather, 
climate and formation of ice in the arctic. Being able to accurately measure 
ocean heat flux is the first step in being able to model and predict these processes 
that influence the habitability of our planet. 

Ocean flux measurement has always been hard to accomplish primarily due 
to the difficulty in sampling the oceans temperature and current structure on 
adequate temporal and spatial scales. Acoustic tomography has a potential for 
overcoming this limitation. In this study, we shall assess the adequacy of a 
tomographic technique for estimating ocean heat flux through simulated arrays 
in an area off the Norwegian Coast. This region is influenced primarily by the flow 
of the Norwegian Atlantic Current. 

The Norwegian Atlantic Current forms a conduit for the poleward flux of 
relatively warm, saline, North Atlantic water along the Norwegian coast shown 
in Figure 1-1. This three Sverdrup! (Sv) flow branches into an easterly flowing 
North Cape Current, and the poleward flowing, West Spitsbergen Current. These 
currents help to mediate the advance and retreat of the ice pack in the Marginal 
Ice Zone. A tomographic array located in this region could in principle allow for 


the monitoring of oceanic heat flux into the Barents Sea and the Arctic basin. 


1 Sverdrup is the unit of mass transport, 1 Sv = 10° meters? sec”! 
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Figure 1-1 Geography: Overview of the Current System of the Nordic Sea 
Region. 


B. APPROACH 

Our assessment begins with a review of some techniques for estimating 
ocean heat flux that have been used in the past. We then assess the tomographic 
method by performing simulation inversions of synthetic acoustic travel time data 


computed from sound speed and current fields derived from climatological 


temperature and salinity information. Separate inversions of the changes and 
differences of reciprocal path travel times will allow us to construct independent 
estimates of temperature and current throughout the volume surrounded by the 
array. From these tomographic products we then construct the zonal and 
meridional estimates of heat flux. We compare the tomographic estimates of 
temperature, current, and heat flux with the actual values of these parameters as 
a means of assessing performance of the technique. The central issue in this 
study is, how well tomography can resolve the temperature and the current fields 
in the ocean, the basic parameters of heat flux calculations. Of particular interest 
is the quality of the current estimate both in magnitude and direction. 

We have selected climatology from the General Digital Environmental Model 
(GDEM) data base with the intent of modeling a region in the upper Norwegian 
Sea off the coast of Norway. The model domain lies in the path of the poleward 
flowing Norwegian Atlantic current. The GDEM data cover two separate 
months that represent seasonal extremes of summer and winter, in a region 
bounded by 5° - 20° East and 67° - 72° North. These two data sets provide the 
ocean environment for our acoustic simulations. 

Our GDEM data base consists of climatological temperature and salinity 
profiles at half degree horizontal resolution. We derive sound velocity and 
potential temperature fields based on empirical expressions using the GDEM 
temperature and salinity fields. We also derive current fields based on geostrophy. 
These derived fields form the truth data in our numerical simulation. 

Additionally, we perform resolution and variance analysis at five depths for 
both the sound speed perturbation and the current field. These statistical analyses 


require no data and they provide us with quantitative measures of how well a 


given array geometry resolves the ocean temperature and current structure. The 
geometries include a large and a small circular array with diameters of 525 and 
250 km respectively. The smaller array contains about 25% of the area of the 
larger array. We would like to determine how much area can be covered by a six 
element array and still adequately measure heat flux. 

We invert the synthetic travel time data obtained from original sound speed 
perturbation and current fields to derive estimates of the sound speed perturbation 
field and current field. We then convert the sound speed perturbation estimate into 
a potential temperature field using a technique explained in the next Chapter. 
Using the potential temperature and current fields, point estimates are made for 
zonal and meridional flux over the model domain. In addition, we compute depth 
integrated estimates and zonal and meridional heat flux through a transect inside 
the array. Estimated values of all the parameters are compared with actual model 

values to assess the quality of each estimate. 

An outline of the remainder of the thesis follows. In the balance of this 
chapter we provide a comparative summary of conventional techniques that have 
been used in ocean heat flux estimates. Additionally, we discuss the proposed 
tomographic method, and the details of how we intend to use it to derive heat flux 
values. Finally, we discuss the environmental factors that are likely to affect the 
performance of a tomographic array in the upper Norwegian Sea. 

In Chapter II, we discuss the forward and inverse problems for both density 
and current tomography. In particular, we discuss how to construct inverse 
solutions, using the Gauss-Markoff estimator, and perform resolution and 
variance analysis that is used to quantify tomographic array performance. We 


also discuss the estimation of heat flux from standard tomographic products, 


sound speed and current maps. In Chapter III, we describe in detail the 
numerical experiment including ocean simulation, array simulation, and the 
creation of synthetic travel time data for the inversion. Chapter IV contains the 
results and discussion of the experiment. Finally in Chapter V, we summarize our 


results and presents our conclusions. 


C. BACKGROUND 
1. Oceanic Heat Flux Measurement 
Heat flux is a quantitative description of the movement of heat through 
a surface from one region to another expressed in units of power (watts). The 
total heat budget for the earth can be viewed as the sum of various terms 
representing sources and sinks. Historically, the approach has been to assume 
that the earth is in thermal equilibrium so that the net energy gain in time is zero. 


This is expressed in Equation (1.1). 


[Net Gain}= [Source] - [Sink] - [Advection ] 
Qner = [Q5- QR) - [Qp + Qp + QH)] - LQya + Qvo J (1.1) 
where: 
Qner= Net heat gain or loss in time. 


Q, = Heat source due to insolation. 


Qe = Heat loss due to reflective albedo. 

Qs = Heat loss due to long wave radiation. 

Q_ = Heat loss due to latent heat of vaporization. 
Q:; = Heat gain due to precipitation. 

Qya = Heat loss due to atmospheric flux. 


Qyo= Heat loss due to ocean flux. 


Average values of the source and sink terms in the heat balance 
equation can be measured with with modest accuracy. Measuring the 
advection terms [ Qy, + Qyo |, is a more difficult problem. Early researchers 
calculated the ocean advection of heat as a small residual of the difference 
between the very large source and sink terms (Sverdrup, 1957). This calculation 
had large uncertainty and could not resolve the controversy over the relative 
importance of the atmosphere and ocean in the flux of heat from the equator. 
The requirement for a steady state and a bounded system did not easily lend itself 
to the study of regional or local heat flux problems. This technique was ultimately 
abandoned in favor of more precise methods. 

The appearance of satellites and the refinement of other in situ sensing 
devices such as current meters led to the development of measurement 
techniques that can be divided into two general classes, namely the direct and 
indirect methods. The technology that made these methods possible marked a 
considerable improvement in the ability to calculate ocean thermal flux. 

a. Direct Methods 

As the name implies, the direct method simply uses direct 
measurement of the components that make up heat flux in the ocean. These 
components are potential temperature @ (K°), and current V (meters sec”'), 
perpendicular to an arbitrary boundary. The general Equation (1.2) is given by 
Bryan (1962). 


Qvo = 


OS mye 


{ pcpo V4 ded (1.2) 
+H 


Heat flux can be broken down into a product of the temperature and velocity 


field perpendicular to a boundary that the heat is being advected through ; where 
L is the length of the boundary, in meters, H is the height in meters, Cy (Joules 
kg 'K°"!) is the specific heat of seawater at constant pressure and p (kg 
meter”) is the density of seawater. 

Since much of the temperature variability of the ocean lies above 
the thermocline, temperature fields can be constructed over broad ocean areas 
based on satellite sea surface temperature maps and several well positioned in 
situ measurements of the water column as shown by Oort and Vonder Haar 
(1976). ‘If satellite data are not available, then direct measurements must be 
made by conventional ocean survey techniques which are laborious, time 
consuming and expensive. 

Velocity fields are more difficult to measure by direct means and a 
variety of techniques have been developed. The most common method is to 
assume that currents are geostrophic. The ageostrophic component is assumed 
to be negligible and for large scale, or climatological studies, this is generally true. 
In local and regional studies, the ageostrophic component may be significant. In 
these cases, other techniques may be necessary if available. Other techniques 
include Eulerian methods involving the deployment of fixed current meters and 
Lagrangian methods that feature sensors moving with the flow that are tracked 
by sonar or aerial reconnaissance. These methods have varying degrees of 
success but each bears a heavy price for coverage and spatial resolution. Thus 
the geostrophic technique remains a popular choice. 

b. Indirect Methods 
Two independent methods have been used for indirect heat flux 


measurement, namely the air-sea exchange and the satellite radiation method. 


The first method is based on the measurement of air-sea fluxes of momentum, 
water vapor, sensible heat and radiation based on bulk aerodynamic and radiation 
equations. The bulk formulas require coefficients that are derived empirically. 
Once the surface fluxes between the air and sea are determined these values are 
integrated over an area poleward of a given latitude. Positive fluxes (air to sea) 
occur at the equator, while negative fluxes (sea to air) occur in the polar regions. 
Heat flux is then calculated by assuming that the deficit in the integrated region 
above a given latitude is exactly balanced by a flux of poleward energy across 
that latitude. A detailed discussion of this method is provided by Bunker (1976). 
An application of this method in the Polar region is given in Mosby (1962). 

A second method is based on the use of satellite radiation data. 
Using this technique, ocean heat flux is calculated as a residual between the net 
incoming and outgoing radiation utilizing the assumptions of Equation (1.1). In 
this method satellite data must be averaged in time and space, and energy 
storage by the ocean and atmosphere is assumed to be negligible. Atmospheric 
heat flux is calculated directly from observations. Ocean flux is then calculated as 
the difference between the total advection and the atmospheric component 
Qvocean = Qvtotat - Qvair- 

The preferred method of computing ocean heat flux is the direct 
method. The method features simplicity and error depends only on how well 
temperature and velocity are measured. The assumptions of a steady state 
system and closed boundaries are unnecessary. Unfortunately temperature and 
in particular velocity are not well resolved over large areas using conventional 
measuring techniques and thus the direct method is costly, time consuming, and 


offers only limited coverage. 


Indirect measurements have the advantage of providing global 
coverage of ocean heat flux. However error can be introduced as in the bulk 
aerodynamic method by the uncertainty of the exchange coefficients that are 
required in these calculations. In addition, many more measurements must be 
taken to compute the various terms required by this technique. This can lead to 
greater error unless the data is carefully averaged. The satellite radiation method 
features excellent areal coverage, but also has some limitations. Diurnal sampling 
and sensor calibration errors, spatial gaps and uncertainties in the value of the 
solar constant can all adversely influence satellite radiation measurements. A 
further drawback is that atmospheric heat flux is not well resolved over the 
oceans due to the sparseness of data, particularly in the southern hemisphere. 

In summary, there are several methods available to compute ocean 
advective heat flux. Resolution varies depending upon the method of choice. 
These measurement techniques fall into two distinct classes, the direct and 
indirect. The direct method is best suited for small scale, high resolution studies. 
The indirect method is most useful in large scale climatological studies. Direct 
methods are accurate but expensive. Indirect methods are more economical but 
have poor temporal resolution and built in bias. 

Comparison of ocean heat flux values computed by these methods 
shows wide variability, especially in lower latitudes. How can this variability be 
reduced? Is there a method that would allow for good spatial and temporal 
resolution without being prohibitively expensive? Tomography may provide the 
answer to these questions. 


2. Tomography Overview 


Prior to detailing tomographic measurement of ocean heat flux, it is 


appropriate to describe some tomography basics. An excellent review of ocean 
acoustic tomography is given in Spindel (1985). 
a. Travel Time 

Eigenrays define fixed acoustic pathways between source and 
receiver. A sound impulse generated at a source travels simultaneously along 
many different eigenpaths to the receiver. Since each path has a different length 
and travels through a different part of the ocean, each eigenray will have a distinct 
and predictable arrival time associated with it. The arrival times of this multipath 


signal are represented in Figure 1-3 for the eigenrays given in Figure 1-4. 
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Figure 1-3 Signal Time Series: Time series of a multipath arrival sequence 
based on the raytrace below. 


As the sound pulses travel they "feel" the environment and are perturbed slightly 
in time by density inhomogeneities and by current along the ray path. By 
measuring the difference between actual arrival time and the arrival time 


predicted for the medium in an unperturbed state, one can infer the sound 
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Figure 1-4 Raytrace: Eigenray trace showing typical multipath pattern from 
the January sound speed profile for the Norwegian Sea region 
from source | to receiver 2. 


speed perturbation field along the path. This is known as density tomography . 

Ocean currents can be inferred by making simultaneous reciprocal 
transmissions between a transceiver pair and examining the difference in time of 
flight along the same eigenpath. Rays traveling with the current will quicken while 
rays traveling opposed to the current will be delayed. This time difference is 
proportional to the the current along the ray path. 

Inverse techniques can be used to reconstruct the four dimensional 
sound speed perturbation and the current field from measured travel time series. 


Further details on inverse theory can be found in a paper by Parker (1977). 
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b. Tomographic Estimate of Heat Flux 
The inversion gives us estimates of the sound speed perturbation 
field and the zonal and meridional components of current. What remains is to 
derive potential temperature estimates so that we can proceed with our estimate 
of heat flux. We do this by assuming that sound speed perturbations are 
proportional to temperature perturbations. This is a valid assumption provided that 
salinity does not vary significantly. We then relate perturbations in sound speed 
to temperature perturbations via the empirical relation giving by Mackenzie 
(1981). Next we take the perturbations in temperature and form potential 
temperature perturbations based on empirical relations given in Unesco (1983). 
Once potential temperature and the components of current are 
found, then it is an easy task to compute point estimates of heat flux. These point 
estimates can then be summed horizontally *rd vertically to yield integrated flux 
estimates. 
3. Environment 
The Norwegian Sea has a number of characteristics that make it a 
challenging location to study heat flux. The waters are complex and diverse both 
physically and dynamically. The poleward flow of warm saline North Atlantic 
water that flows along the Norwegian coast is a sharp contrast to the cooler, less 
saline waters of the Arctic Basin. In this environment, meanders, jets rings and 
gyres are not uncommon. Additionally, the ocean's temperature is also quite 
variable. Sharp thermal gradients are formed at the confluence of the different 
water masses. These features will provide a good test of the tomographic method. 
The location of the region where our numerical experiment is conducted is 


shown in Figure 1-5. 
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The environment plays a major role in the way heat is transported in the Arctic 
basin, and it also critically affects acoustic propagation in the region. Seasonal 
changes modify the water mass structure and directly influence acoustic 


propagation creating two quite diverse regimes. Bathymetric features 
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Figure 1-5 Simulation Location: Location of the Tomographic Array for the 
measurement of Heat Flux. 


Seasonal ice formation and breakup processes add complexity to the heat budget 
such as the shelf break have a controlling influence on the mean flows and the 
acoustics of the region. Fronts are created at the boundaries of dissimilar water 
masses. Each of these factors will have an impact on the deployment of a 
tomographic system. 

a. Dominant Flow 

The Norwegian Atlantic Current and the Norwegian Coastal 
Current constitute the dominant flows in the region. Both of these currents are 
northeasterly, hugging the Norwegian coastline. The Norwegian Atlantic 
Current is an extension of the Gulf Stream that makes its way through the 
Faeroe-Shetland Channel importing relatively warm saline water from the North 
Atlantic. The temperature of this water drops from 8° to 3° Celsius in the winter 
and 10° to 6° Celsius in the summer by the time it reaches the Spitsbergen area. 
The Norwegian Atlantic Current is roughly 80 km wide, 500 m deep, flows at 20- 
40 cm sec’! and is responsible for a transport of approximately 2-3 Sv. This 
current is extensively reviewed by Saelen (1962). 

The No.wegian Coastal current consists of a wedge of slightly less 
saline water representing effluent from the Baltic and runoff from the myriad of 
fiords that are scattered along the Norwegian coast. Flow of the Norwegian 
Coastal current is considerably less, and has been neglected in previous heat 
budget calculations. 

b. Bathymetry 

These flows are heavily influenced by the bathymetry of the region, 

details of which are described in Hurdle (1986). As the Norwegian Atlantic 


Current flows north it traces a path directly above the sharply sloping shelf break 


where it becomes the northerly flowing branch of the West Spitsbergen Current. 
The West Spitsbergen Current seems to elude topographic influence from time to 
time and has been observed to branch into westerly and southerly flows. The 
easterly flowing North Cape Current is the exception crossing the shelf break and 
flowing into the shallow Barents Sea basin. 

Although under the influence of bathymetry, the Norwegian 
Atlantic Current is by no means a rigid flow. This current can meander and spin 
off mesoscale eddies. These features can be advected northward and play a role 
in the melting of the ice sheets in the marginal ice zone. A description of these 
eddies created by this flow is given in Saelen (1963). 

c. Fronts 

Fronts can be expected to exist at the confluence of different water 
masses and these features complicate sound propagation. Fronts vary in strength 
and thus their acoustic effect. In this region, they are permanent features that 
vary in intensity on a seasonal basis with maximum strength in winter. Two 
fronts lie in close proximity to this study. The Polar Ocean Front is the strongest 
extending from Spitsbergen southward to Bear Island where the gradient is most 
intense at 3-4° Celsius per 15 km. The front then turns eastward into the Barents 
Sea and forms a line of demarcation between the icebound regions to the north 
and the ice free waters to the south. A weaker extension of the Polar Ocean Front 
branches westward from Bear Island into the Norwegian Sea. Acoustic paths 
can be quite variable across fronts making raypath identification difficult. 

d. Ice 
Ice is a physical barrier that can alter the acoustic path and the 


travel time and thus has an impact on tomographic system performance. Two 


different regimes are apparent. The first involves sheet ice which alters the 
acoustic reflection characteristics at the surface boundary and the second is free 
floating ice that may temporarily affect the acoustic path. 

The advance of sheet ice in the Barents and Norwegian Sea is 
mediated by the flow of warm water brought northward by the Norwegian 
Atlantic Current. The confluence of warm and cold water in the region forms a 
frontal boundary called the Polar Ocean Front which marks the southern 
boundary of the ice in this region We have located the tomographic array well 
south of this front, thus avoiding the effects of sheet ice altogether. 

Concern over free floating ice is minimal due to the northerly flow 
of warm saline North Atlantic water along the Norwegian coast. Should large ice 
flows be carried out of the Barents by the East Spitzbergen current, a potential 
disruption in acoustic paths could occur. These cases would degrade system 
performance if acoustic pathways were obstructed however, the effect would 
only be transient in nature and disappear as the ice is advected away. 

e. Acoustic Properties 

Seasonal changes in the region bring about quite remarkable 
changes in acoustic propagation in the region. The mean background sound 
speed profiles are shown in Figure 1-6 for the two monthly climatologies 


selected for study. 
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Figure 1-6 Acoustic Environment: Plot of Background Sound Velocity 
Profiles for January and August 1989. Error bars represent one 
standard deviation. 
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Both climatologies show almost identical deep sound channels located at 
approximately 1000 meters. The near surface profile are very different however. 
The January 1989 profile shows a local sound speed maximum at roughly 400 
meters which establishes a surface duct. This ducting is created by the cooler 
surface waters that characterize the winter months. The August 1989 profile in 
contrast shows a decreasing sound speed profile from the surface to the deep 
sound channel axis, creating a very strong channeling environment in the first 
kilometer. This channeling environment is a result of the warmer surface waters 
that characterize the summer months. These two contrasting acoustic regimes 
will dictate the depth of the array transceiver and the type of eigenrays available 
for use in this study. 
f Summary 

In summary, the environment is expected to have an important 
impact on tomographic system performance, and does present some serious 
system complications. Most of these effects are difficult to model and are treated 


as limitations to scope in this initial study. 
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Il. PHYSICS OVERVIEW 


In this chapter we present an overview of the “forward problem” associated 
with density and current tomography. The “forward problem” specifies the 
relationship between the observables, (i.e. travel times), and the unknown sound 
speed perturbation and current fields, we are estimating. We examine the use of 
a stochastic inverse method that allows us to form estimates of these unknown 
fields from the observed data. In addition, we discuss statistical measures of 
system performance. These measures include resolution length and variance. A 


derivation of heat flux from current and sound speed estimates is also included. 


A. TOMOGRAPHIC FORWARD MODELING 
1. Time of Flight 

In our acoustic model, we require an accurate representation of the path 
that a sound wave front takes as it travels through the ocean. The path is 
important because by predetermining the path we can calculate expected travel 
time of the signals between transceivers in an array. A convenient approach is to 
model sound propagation using ray theory. We employ ray theory for two basic 
reasons. Aside from it's simplicity, the ray approach offers a good physical 
representation of the multipath propagation of sound in an acoustic duct. The 
theory breaks down near caustics and at turning points, consequently corrections 
must be made. These limitations are not significant enough to lessen the 


usefulness of this approach. The time of flight of a sound pulse along a given 


eigenray connecting a source and receiver can be evaluated as follows: 


ds 


T,+6T, = i es oe (2.1) 
Cotic + Ve S 


path; 


a. Density Tomography 
Equation (2.1) reveals that the time of flight along any given raypath 
is inversely proportional to the sound speed and current integrated along the ray 
path. Upon linearization, the change in travel time resulting from the presence of 


sound speed perturbations and the flow field along the path and can be expressed 
as 
Se) | (2.2) 


path; 


Since l&c |>>1 Ve § | we can approximate Equation (2.2 ) by 


oT. = - bc ds (2.3) 


3 
Co 


path; 


It can be seen that travel time perturbation is proportional to the sound speed 


perturbation summed along the ray path. 
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b. Current Tomography 
Current is related to the difference between travel times along 


reciprocal paths. A subtraction of reciprocal travel times yield 


8T, = +t. -—_ +. s ( 2.4) 
Cotdc+VeS Cy+dc- VeS 
path; 


upon linearization we get 


(2.5) 


path; 


2. Experimental Noise 

The forward model is not complete without inclusion of experimental 
noise. Travel time measurements are contaminated by noise that contains two 
components, namely measurement noise and model error. Measurement noise 
results from the lateral] displacement of the transceiver element caused by 
internal waves or tidal effects. Model noise is due to the limits on numerical 
accuracy and the assumptions that led to the linearization, and discretization of 
the forward problem. In practice the effects of measurement noise can be 
circumvented by currecting for transceiver motion or employing various signal 


processing techniques. A technique that is commonly used is to average the 


travel time measurements over several transmission cycles. These techniques, as 


well as others, are discussed in Spindel (1985). We represent these effects by 


using a random noise generator to add a noise component to our travel time data. 


B. A DISCRETE STOCHASTIC INVERSE METHOD 
1. Discretization of an Ocean Volume 

In our solution to the inverse problem, we have chosen to discretize the 
simulated ocean region by dividing it into 605 equal volume boxes measuring 50 x 
50 x 0.2 km each. Within each box the sound speed perturbation and current are 
assumed to be constant. This discretization results in a total of 11x 11 squares 
in the horizontal and 5 layers in the vertical. This allows us to cast the continuous 
integrals in (2.3) and (2.5) as discrete sums. Placing all the data in column 


vectors, we can re-express the forward problem in matrix-vector form as 


6T =A dc+e (2.6) 


dT = B V+ e' (2.7) 


where 6c and V are vector parameterizations of the sound speed and current 
fields respectively. A and B are linear operators that act as transfer functions 
relating the unknowns to the data. The vectors e and e' represent the 
experimental noise. We assume that these errors are uncorrelated with the signal 
and have a known variance and a zero mean. 
2. Gauss-Markoff Estimate 

Equations (2.6) and (2.7) represent a highly underdetermined system, 
that is there are many more unknowns than equations or data. There is infinite 
number of solutions satisfying these equations. The best solution however, should 


fit the data as well as the statistical a priori information. A useful method that 
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requires the optimal solution to have minimum mean square error is given by the 
Gauss-Markoff theorem. This method states that the optimal estimate can be 


A 
evaluated, assuming &c and V are uncorrelated, as 


5¢- C. ATC, OT (2.8) 


V = GB'Ce 8T- (2.9) 


where 


Ce = Cy,-(C,A™ (AC, AT+ Cey(C,, A) 2.10) 


& dc 4+ 
C,= Cy-(Cy B™ BC, B™+ Ce)’ (CyB (2.11) 


are the error covariances of 5¢ and Vv respectively, and Ce and Ce’ are the 
covariance matrices of experimental noise ¢ and e' respectively. An appealing 
property of the Gauss-Markoff estimator is that it does not require the shape of 
the probability density functions of the unknowns. 
3. Covariance Functions 

In order to reconstruct the sound speed perturbation field and current 
field using the Gauss-Markoff method we need to specify the covariances of the 
fields. It is through the covariance matrices that we are able to specify “a priori 
knowledge” concerning the unknowns. In the inversion of the data, we assume 
that the sound speed perturbation field and the zonal and meridional current field 


are uncorrelated with each other. Furthermore we assume that all these fields are 
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Statistically homogeneous and their correlation functions have gaussian shapes, 


that is 
Cc _ x2 -| (AX)2 2 a 
OV,,(Ax,Ay,Az) = 05, exp | (F-) + (72 (2.12) 
Cov, (Ax,Ay,Az) = Oy (exp” ee ) (oy (7 5p (2.13) 
Covy (Ax.dy.de) = Oy (exp [2 )24(F cL) + (72 ai (2.14) 


The covariance function is inexact, but a good approximation if the 
correlation lengths L,, Ly, and L, approach the true correlation lengths of the 
ocean environment. that correspond to the true ocean environment. The 
Ax,Ay,Az values represent the separation between two points in the given 


direction. The quantities 05,,0y, and oy y are the rms values of the respective 


unknowns. Values for 65, are roughly 5 meters sec! while oy. values are around 


0.5 meters sec-!. The larger the correlation length, the more information we 
specify to the estimator. In the inversions we used horizontal lengths of 50-100 
km and vertical length of .2-.6 km. These values lie at the upper limit of size for 
mesoscale features in the region. 
4. Error 
Error in the estimate, denoted by subscripts € and N in Equations (2.8) 
and (2.9), can be broken down into two Statistically independent components. 


These components are the bias, 
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b=< 8 > - Sc (2.15) 


dse Vs 2V (2.16) 
and random error 
A( 8¢)= 8&-< &> (2.17) 
AV =V-<V> (2.18) 


where < 6c > and< Vv > are the expected value of Sc and Vv respectfully. 

Bias is attributed to systematic error, which is a function of sample size. 
As the size of the sample data increases, bias is expected to decrease and vice 
versa. Random error results from the randomness of the data itself, which 


contains measurement error and model error. The total mean square error is 


= ‘ T 
C- Ce +<bb' > (2.19) 
C.- Cas + <dd™> (2.20) 
The covariances of the randomness are: 


-1 
Cis =C A'CeAC. (2.21) 


Cys = C,B'CeBC, (2.22) 
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5. Resolution 
Having derived the statistical measures of error we now derive 
measures of system resolution. The following discussion of system resolution 
starts with an application of a singular value decomposition as described by Chiu 


(1987). The operators in Equation (2.21) and (2.22) can be decomposed as 


[co"A CY] = UAV’ (2.23) 
[ce"B C?] = QPS™ (2.24) 


where the diagonal elements A, of matrix A and $, of matrix P are the 


s. of the matrices U , V 


associated singular values, and the columns u., v. and q;, S, 


and Q , S represent the left and right singular vectors defining the parameter and 
data spaces. 
Following the derivation given by Chiu et al (1987), the minimum 


mean square error estimate can be expressed in terms of the singular vectors as 


“1/2 eA Mi Tpw-l/2 
Cy dC = a [us C. oT | V; (2.25) 
o + 


i=l 


The solution for each estimate is obtained by taking a weighted sum of the base 


vectors V; and §;. The diagonal elements A , and >. are analogous to signal to 


noise ratios. For example, when an? or o; >> 1, the signal dominates, and in the 
opposite case, noise dominates. The linear optimal estimator acts to downweight - 
the noise while passing the signal. 


Following the discussion of Chiu (1987), we define a symmetric 


matrix R as the resolution matrix such that 


Ry= VAC - A’)! AV" (2.27) 


Ry = SP(1-P?)'PS" (2.28) 


The rows and columns in this matrix correspond to the resolution kernel, which 
represents the best estimate of a delta function. Substituting the resolution matrix 


into Equations (2.27) and (2.28) and performing some manipulation, we come up 


with 
Cre Cx (1- Ry Cx (2.29) 
Cn- Cy’ (I- Ry )Cv" (2.30) 


These Equations shows that for the minimum mean square error 


estimate, the relationship between error and resolution is linear. As resolution 


increases, error decreases and vice-versa. When R=I, Ce= 0, the ideal 
resolution situation exists. 

As a means of describing the system resolution, we adopt Chiu's 
(1987) measure of horizontal and vertical resolution length H and V respectively. 
For each unknown field, i.e. sound speed perturbation or current, these quantities 
represent the horizontal and vertical distances for which approxiniately 70% 
(+1 © ) of the total energy of the resolution kernel associated with a given 
location is confined. The minimum resolution lengths in the three spatial 


coordinates is calculated by 


H, (ix, iy, iz) =2 oY Y [Gx-ixddx]? ri (jx.iy.jz) E;’ (2.31) 

JX JY JZ 

ee nxMynz, 27. ss : 

Hi fix,iy,iz) =2 by oy [Gy-iy)dy }* ri Gx. jy.jz) E} (2.32) 

JX JY Jz 
Vix,iyiz) = 24 | SY ¥ [Gz-iz)dz]* 1 (jx,y.iz) E;' (233) 

JX JY Jz 

where 
, mxnynz 2, _ 

E;' = 2 z z ri (jx, jy,jz) (2.34) 
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where jx, jy, jz are three dimensional box indices and 4g is the resolution kernel. 


H xX» H y and V represent the minimum resolution length in each spatial 
coordinate. An important point to note is that RMS error and resolution length do 
not depend on the environmental data 5T, but rather on the sound speed 
perturbation covariance matrix Coc or the covariance matrix of the current Cy 
and the array geometry. 
6. Heat flux estimate and error 

We can estimate sensible temperature from the sound speed field by 
using an empirical relation given in Mackenzie (1983). Sensible temperature can 
be related to potential temperature by an empirical relation given in UNESCO 
(1981). We use these relationships to compute the potential temperature from 
the sound speed field throughout the model. Heat flux estimates can be 


constructed from @ and V as follows 
F= pc, 6 %, (2.35) 
F pc, 6 yy (2.36) 


where F, PF, and V,, v, are the zonal and meridional components of heat flux 
and current respectively. 

The error in the heat flux estimate is related to the error in the individual 
parameters. Assuming the individual errors in 6 ; v. , and v, are uncorrelated, 


we can determine the variance of each estimate by letting 


€,=0-8 (2.37) 


ane = V, - Vx (2.38) 


co = V, - Vx (2.39) 


be the individual errors. The mean square error of each flux estimate is 


ce, > =< (f, -fy?> =p Cp<0V,-8Vx> (2.40) 


<e, 2, =< (f, - fy)?> =p Cp< O Vy - 6 Vy> (2.41) 
y 


Substituting the error represented in (2.37) , (2.38) , and (2.39) into Equation 


(2.40) and (2.41) and collecting the lowest order terms we have 


A 
o 
Vv 

I 


= [ 0? <ey *> +Vx,2<e,"> | (2.42) 


A 
o 
pens 
Vv 
WW 


[ 6? <ey *> + Vy? <e,"> | (2.43) 


Assuming that the errors are small one can replace the actual values 0, V,, and 


V, in (2.42) and (2.43) with the estimates 6,9, ,and_ ,. 


C. SUMMARY 

In summary, we have presented an overview of the forward problem of 
density and current tomography. We have shown the relationship between the 
travel time data and the sound speed perturbation and current fields we are 
estimating. The derivation of the heat flux estimate and its mean square error 


were presented. We then reviewed the Gauss-Markoff technique which enables 
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us to invert the data for optimal estimates of the unknown fields. We discuss 
some useful statistical measures of system performance.These are the horizontal 


and vertical resolution lengths and mean square error. 
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III. NUMERICAL OCEAN AND ARRAY SIMULATIONS 


A. SYNTHETIC OCEAN 

The synthetic ocean represents the "truth" in this numerical experiment. A 
realistic Norwegian Sea environment is simulated through the use of GDEM 
climatology. 

1. Spatial Resolution 

Selection of spatial resolution becomes a tradeoff between computer 
resources available and the desire to express the fine spatial detail to adequately 
describe the flow in the region. The simulation experiment was performed on a 
VAX GPX minicomputer. A reasonable tradeoff between resolution and 
computation time limited the model domain to a region of 550 x 550 x 1 km?. The 
model domain was subdivided into 605 equal-volume boxes each measuring 50 x 
50 x .2 km. This arrangement allows for simulation of a 302,500 km? ocean 
volume at a resolution high enough to include most mesoscale variability in the 
upper | km. 

Vertical resolution of the synthetic ocean is an important consideration 
given the baroclinic nature of the temperature field and the moderate vertical 
shear seen in the geostrophic current data. The strongest vertical temperature 
and current gradients are located in the near surface regime. A 200 meter depth 
increment was chosen to provide the optimal balance between adequately 


defining the field and keeping the problem tractable on a minicomputer. 
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2. Bathymetry 
Due to the unavailability of precise bathymetry and the uncertainties of 
modeling absorption and reflection with the bottom in the region, we have 
restricted our study to only non-bottom interacting rays. 
3. GDEM Data Base 
The model ocean is dervived from GDEM climatology made available 
by Fleet Numerical Oceanography Center. GDEM is a synthetic data base of 
temperature and salinity profiles that cover 30 depths irregularly spaced from the 
surface to 5000 meters. Synthetic profiles of temperature and salinity are based 
on regression equations built from a data base of over four million observations. 
GDEM climatologies are available in monthly, seasonal or annual time frames. In 
this experiment, we choose two mean climatologies for the months of January 89, 
and August 89, representing annual extremes in temperature. The irregularly 
spaced GDEM data given in geodetic coordinates were mapped into a regularly 
spaced Cartesian system in the model domain using cubic splines. The 
climatological data base contains information only on the average structure of the 
ocean. Therefore, some mesoscale and sub-mesoscale variability that are present 
in the real world are not represented in the synthesized ocean, 
4. Model Parameters 
Once temperature and salinity data were obtained on a regular grid, 
other model parameters were computed at each grid point. Sound speed was 
computed using the empirical formula given by Mackenzie (1981) and then 
separated into a background, Co(z), and a perturbation field 5c(x,y,z). Potential 
temperature was calculated using the empirical formula given in UNESCO 


(1983) using the GDEM temperature and salinity data. The actual sound speed 
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perturbation field and potential temperature fields for 5 levels of the model for 
January and August 1989 are presented in Figures 3-1 and 3-2. 

Geostrophic current was also computed using the GDEM temperature 
and salinity data. The geostrophic approximation is used to derive the velocity 
field because no current fields were available, in the data base and geostrophy 
represents a physical way to generate synthetic velocity fields. Figure 3-3 shows 


the derived current fields for January and August 1989 respectively. 
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Figure 3-1 Sound Speed Perturbation Field: Contouring of the actual 
January and August 89 sound speed perturbation fields at three 
depth levels. Units are in meters sec”! 
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August 89 potential temperature fields in degrees C at three 


depth levels. 
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Figure 3-3 Current Fields: The “actual” January and August 89 geostrophic 
velocity field at three depth levels. The magnitude of the current 


ranges from .1-.5 meters sec’ at these depths. 
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5. Acoustic Environment 

Both climatologies show a deep sound channel located at approximately 
1000 meters (see Figure 1-7). The near surface profiles are very different 
however. The January 1989 profile shows a sound speed maximum at roughly 
400 meters which establishes a surface duct. This ducting is created by the 
cooler surface waters that characterize the winter months. The August 1989 
profile in contrast shows decreasing sound speed from the surface to the deep 
sound channel axis, creating a channeling environment. This channeling 
environment is a result of the warmer surface temperatures that characterize the 
summer months. These two contrasting acoustic regimes will govern the optimal 
depth of the array sources and receivers and the type of eigenrays available for 


use. 


B. ARRAY GEOMETRY 

Two six, element arrays are simulated in this experiment, and they are 
depicted in Figure 3-4. The large array has been given the name Steven and the 
small array, Harry. These geometries were selected to answer two basic 
questions. The first one is, how much area can be covered by a six element array 
whose purpose is the measurement of heat flux? Another question is how does 
system resolution vary with array size. 

These two array geometries should in principle have advantages under 
different circumstances. Steven should be better suited to cover a larger area 
with less resolution, while Harry would be better suited to cover a smaller area at 
greater resolution . Both geometries are compared in a quantitative sense and 


many more arrangements are possible. 
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Figure 3-4 Array Geometry: Two 6 element arrays, Steven (large) and 
Harry (small) , simulated in this experiment. 


C. RAY TRACING 

With the background sound speed profile determined array geometries 
selected ray tracing can proceed. Eigenrays are calculated using a raytrace 
algorithm based on a 4th order Runge-Kutta numerical integration technique for 
selected horizontal configurations and transceiver depths. In this experiment we 
identified all eigenrays that had a launch angle lying between -15° and +15°. 
Eigenrays having unique trajectories were selected to comprise the data set for 
inversion. This provided between 3 and 10 eigenrays for each transceiver pair. 
Sample eigenray traces are shown in Figures (3-5) and (3-6) for January and 


August 89 climatology. 
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Figure 3-5 Raytrace: Eigenray trace for a transceiver pair from array Steven 
from January 1989 climatology. 


OETA tem 
eee 


ten ees ¢€ BW 40 C8 fe 


ewe eT 


18@ 1280 «6148 «6308 «6188 «288 «228 248 288 28: 
RANGE (Ke) 


Figure 3-6 Raytrace: Eigenray trace for a transceiver pair from the Steven 
array from August 1989 climatology. 
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D. SYNTHETIC TRAVEL TIME DATA 
Synthetic travel time data are computed for the two array geometries using 
the actual sound speed perturbation and current fields We add to these signals a 


random error component to simulate experimental noise. 


E. INVERSION 

Our intention is to assess the tomographic method on it's own merits Thus in 
our inversion we specify no physical or dynamical constraints. Travel time data 
are inverted using two specialized software programs, one for density and the 
other for current. The output of these programs are the estimate of sound speed 
perturbation and the zonal and meridional components of current. These 
estimated fields are compared to the actual fields at each level in the model. In 
addition , the reconstructed perturbation fields are added to the background field 
to obtain the total sound speed field. Estimates of potential temperature are made 
using the total sound speed field as previously described. Point estimates of zonal 


and meridional heat flux are then computed. 


F. RESOLUTION 

The performance of an array can be assessed using statistical means 
independent of any data. A system resolution analysis can help to determine the 
appropriate geometry for heat flux measurement. We can calculate the horizontal 
and vertical resolution lengths of each array using the method described in the 
previous chapter. 

A desirable array is one which can adequately resolve both the temperature 
and current fields, and yet can cover as much area as possible. How much area 


can a six element array cover and still adequately resolve temperature and 
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current structure in the Norwegian Sea? Which measurement, temperature or 


current will have the most stringent resolution requirements. for the determination 
of heat flux Answer these questions are given in a resolution analysis presented 


in the next chapter. 
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IV. RESULTS AND DISCUSSION 


In this chapter we discuss the results of our investigation. First we show a 
resolution analysis which quantifies tomographic array performance. Next, we 
compare the actual error in the estimated fields to the root mean square error 
predicted by the Gauss-Markoff estimator. A direct comparison of the actual 
and estimated fields is then followed. We do this for potential temperature, 
current, and then heat flux. We consider our assessment to be a conservative 
one since we have used rather inexact a priori information and high noise levels 
for the monthly estimates. 

The two climatologies used in this simulation provide a formidable test of the 
tomographic system. The temperature fields for each month contain some rather 
fine horizontal structure. For example, both climatologies contain a warm core 
eddy feature, that is expressed at 500 meters and below. At depths above 500 
meters a frontal feature is present in the North-West quadrant of the model 
domain. The current field is also complex. There is a anticyclonic gyre 
superimposed on a warm core feature and rather sharp oscillations are apparent 
in the both zonal and meridional components of the current field. The meridional 
component of the current is dominant while the zonal component is weaker. The 
question of how well these fields can be estimated depends on how well 


tomography samples the ocean volume. 
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A. RESOLUTION 
1. Sound speed perturbation field 

Resolution length defines the smallest feature that can be resolved by an 
array. Horizontal resolution depends on array geometry. Vertical resolution 
depends on the characteristics of the acoustic environment. Figure 4-1 through 
4-3 show the minimum horizontal and vertical resolution length for the January 89 
sound speed perturbation field. These results are summarized in Table 4.1. 
Horizontal resolution lengths for the sound speed perturbation field averages 
about 120 km for array Steven while for array Harry it averages 30 km. The 
improved resolution for the Harry is due to the compact design. Horizontal 
resolution decreases slightly with depth for both arrays. Harry, by virtue of 
compactness provides better resolution overall and seems to do slightly better in 
the meridional direction. Vertical resolution is about the same for both 
geometries, with a minimum resolution length of 440-480 meters. Vertical 


resolution depends heavily on the trajectory of the eigenrays . 


TABLE 4.1 Summary of the resolution analysis of the January dc field. 


Resolution Length Array Steven Array Harry 
Zonal (km) 120 30 
Meridional (km) 120 25 
Vertical (m) 400 240 


A look at horizontal resolution for the August case reveals that the 
horizontal resolution is about the same; however, vertical resolution is 30-40% 


poorer. The poorer vertical resolution is due to the channeling nature of the 
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acoustic environment. A way to improve the vertical resolution of the August 
case would be to couple tomography with satellite measurement of sea surface 
temperature (SST). The inclusion of satellite data is suggested as a topic for 
future study. 

In summary the horizontal resolution array Harry is four times better than 
array Steven for the sound speed perturbation field. The vertical resolution was 
roughly 15% better for Harry when compared to array Steven. Horizontal 
resolution length is roughly 20% of the array diameter for array Steven while 


12% of array diameter for array Harry. 
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Figure 4-1 Resolution: Minimum Zonal Resolution length (km) of the January 
89 sound speed perturbation field for two array geometries. 
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Figure 4-2 Resolution: Minimum meridional resolution length (km) of January 


89 sound speed perturbation field for two array geometries. 
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Figure 4-3 Resolution: Minimum Vertical Resolution length (m) of January 
89 sound speed perturbation field for two array geometries . 
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2. Current field 


Resolution of the January 89 current field is presented for the two arrays 


in Figures 4-4 through 4-9. A summary of these results is given in Table 4.2. In 


comparing the resolution of the zonal component of current for both arrays in the 


zonal direction, we observe in Figure 4-4 that Harry does significantly better with 


a resolution length of 35 km. Zonal resolution in the meridional direction is 25 km 


for Harry and remains at 125 km for Steven as is shown in Figure 4-5. 


TABLE 4-2 Summary of Resolution Analysis for current. 


Resolution Length 
Zonal Current 

Zonal Direction km 
Zonal Current 
Meridiona] Direction km 
Zonal Current 

Vertical Direction m 
Meridional Current 
Zonal Direction km 
Meridional Current 
Meridional Direction km 
Meridional Current 


Vertical Direction m 
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Figure 4-4 Resolution: Minimum resolution length (km) of the zonal 


component of January 89 current in the zonal direction for two 


arrays . 
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Figure 4-5 Resolution: Minimum resolution length (km) of the zonal 
component of current in the meridional direction for the two arrays 
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Figure 4-6 Resolution: Minimum resolution length (km) for the meridional 
component of current in the zonal direction for two arrays . 
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Figure 4-7 Resolution: Minimum resolution length (km) for the meridional 
component of current in the meridional direction for two arrays. 
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Figure 4-8 Resolution: Minimum Vertical resolution length (m) for the zonal 
component of current for two arrays. 
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Figure 4-9 Resolution: Minimum Vertical resolution length (m) for the 
meridional component of current for two arrays. 


Array Harry is by far the superior array when it comes to resolution. The 
small array outperforms it's larger competitor by a 2 to 1 margin in the horizontal 
and 16% in the vertical. It is obvious that increased resolution is at the expense of 
coverage. The compact array covers only about 25% of the area of the larger 
array. 

Figure 4-10 is a plot of the tradeoff between minimum resolution length and 
array diameter, plotted for both sound speed perturbation and current. Apparent 
from this plot is that the sound speed perturbation field and the current field have 
different minimum resolution lengths. This is most pronounced near the array 
center. The current field has a larger minimum resolution length indicating that the 
array has more difficulty resolving these fields. The minimum resolution length for 
sound speed and current is roughly 20% of the array diameter for array Steven. 
There is a difference in resolution of the two unknowns in the smaller array. The 
minimum resolution length of sound speed is 25 km while for current it is 35 km. 
This means that for array Harry the minimum resolution length is 10% of array 
diameter for sound speed and 14% for current. Although the current field is more 


difficult to resolve we can improve resolution by reducing array size. 
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Figure 4-10 Resolution: Tradeoff between array diameter and resolution for 
sound speed perturbation and current fields. 


Resolution length tells us whether an array has an adequate resolution to 
estimate the fields. An appropriate array is one that has a resolution length that is 
equal to or smaller than the scale of spatial variability one wants to resolve. Since 
our data field has a horizontal scale of 50 km, we would require a minimum 


resolution length of at least 50 km to adequately resolve the unknown fields. 


Array Harry with its small minimum resolution length satisfies this requirement. 


If the resolution length is greater than 50 km, the estimated field will be a 
smoothed approximation of the actual field. The estimator will act as a low pass 
filter. This will become clearer as we examine our error analysis results in the 


next section. 


B. ERROR 
Error in the estimate is viewed in two different ways. First we examine 
actual errors by comparing the actual vs. estimated fields for temperature current 
and heat flux. We then compare the actual error with the predicted root mean 
square error of the estimated sound speed perturbation and current fields based 
on the Gauss-Markoff estimator. 
1. Actual vs. Estimated Fields 
a. Theta 
Actual error in the estimate of theta was smallest. The actual and 
- estimated fields from array Steven are compared in Figure 4-11. Comparison of 
the fields shows that tomography does a good job in resolving the temperature 
structure both horizontally and vertically. The estimate smooths some of the fine 
structure present in the temperature field, however this effect is subtle. Actual 
error is less than 1° C inside the array. Array Harry demonstrated similar 
performance. 
In the August 89 case horizontal resolution was adequate but 
vertical resoluticn was poor. This is due to the lack of surface turning rays 


present in the acoustic environment. 
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Figure 4-11 Actual vs Estimated Theta Comparison of actual vs. estimated 


January 89 theta fields. 
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Potential temperature estimates based on the compact geometry 


showed minor improvements. Estimates of the summer potential temperature 
field were not as good as the winter estimates. This again is attributed to a 
poorer vertical resolution. 

b. Current estimates 

Current is the parameter having two degrees of freedom, 
magnitude and direction, and a higher spatial variability and thus it is expected to 
be more difficult to estimate. Figures 4-12 and 4-13 show the actual vs. 
estimated fields for January 89 current for both arrays. The estimated fields are a 
smoothed version of the actual flow fields. Array Steven has the largest errors in 
magnitude and direction while array Harry resolves the flow quite well. 

The large error in Steven can be explained in part by the nature of 
the a priori information we have supplied the estimator. Our estimator requires 
thai we specify the covariance of the parameter field. We have assumed that 
each field is independent, isotropic, and has moderate spatial correlation. The dc 
parameter meets these criteria, however the current fields are quite different. 
Not : that in Figure 4-12 the actual flow field is highly anisotropic and has a small 
spaual correlation. Steven does a very poor job in estimating the current 
espt cially the meridional field due to the inaccuracy of the a priori information we 
have specified in addition to having inadequate resolution. Both arrays tends to 
under-estimate the meridional field while over-estimating the zonal field. This is 
due to a leakage or redistribution of energy from the meridional to the zonal field 


caused by the inexact assumption of isotrophy. 
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One way to overcome an inaccurate specification of a priori 
information is to reduce the array diameter, and this improves the horizontal 
resolution. Despite the inaccuracy of the a priori information array Harry does a 
good job in resolving the current . 

The current field has a higher variability. The smoothing seen in the 
current estimated by array Steven is due to variations in current on a scale 
smaller than the minimum resolution length for that array. The smoothing of the 
current is significantly reduced for array Harry because the resolution length 
approaches the scale of the current variation. Note that temperature is resolved 
well in both arrays. The implication is that we must design an array to fully 


resolve the current structure as a minimum requirement. 
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Figure 4-11 Actual vs Estimated Current: Comparison of the Actual vs 
estimated current for array Steven. 
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Figure 4-12 Actual vs Estimated Current 


estimated current for array Harry. 
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2. Predicted RMS Error 
a. Predicted error of the sound speed field. 
A comparison of actual and predicted error is presented in Figures 
4-13 and 4-14 for both arrays. The estimator predicts less than a 50% error in the 
sound speed perturbation field for array Steven and Jess than 25% error for array 
Harry. Actual error for both arrays was less than predicted indicating that the 
actual sound speed perturbation field has a much longer correlation length than 
that specified. 
b. Predicted error of the current estimate 
The predicted error for the current field is presented in Figures 4-15 
through 4-18. The estimator predicts an error of less than 50% for array Harry 
and less than 75% for array Steven. The larger error in Steven is associated with 


with a poorer resolution. 
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Figure 4-13 Sound Speed Error: Comparison of actual vs predicted error in 
the estimate of sound speed perturbation in percent by array 
Steven 
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Figure 4-14 Sound Speed Error: Comparison of actual vs predicted error in 
the estimate of sound speed perturbation in percent by array Harry 
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Figure 4-15 Current Error: Comparison of actual vs predicted error in the 
estimate of zonal current field in percent by array Steven. 
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Figure 4-16 Current Error: Comparison of actual vs predicted error in the 
estimate of meridional current field in percent by array Steven. 
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Figure 4-17 Current Error: Comparison of actual vs predicted error in the 
estimate of zonal current field in percent by array Harry. 
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Figure 4-18 Current Error: Comparison of actual vs predicted error in the 


estimate of meridional current field in percent by array Harry. 
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C. HEAT FLUX ESTIMATES 

Now we examine heat flux estimates within the array from three different 
perspectives, namely point estimates, depth integrated estimates and heat flux 
through a transect within each array. 

1. Point estimates. 

In Figure 4-19, we compare the actual and estimated January 89 zonal 
flux values at the top three layers, 100, 300, and SOO meters respectively. The 
estimates are based on arrays Steven and Harry, respectively. Comparing the 
fields, we note that array Steven tends to smooth the actual fields. In addition it 
also tends to under-estimate in the top layers while over-estimating the deeper 
layers. There is a small area near the center of the array where an over- 
estimation occurs. This is due to the inability of Steven to resolve the fine spatial 
scale. 

Zonal flux estimates for array Harry compare more favorably with the 
actual values. Although this array covers a much smaller area the estimates are 
very precise, resolving well the direction and amplitude of the flux fields as shown 
in Figure 4-20. 

Meridional estimates are compared to the actual flux values in Figure 
4-21. Considerably more structure is evident in the meridional fields with a wider 
range between maxima and minimum. Most of the complex structure is 
concentrated above 500 m and tends to mirror in character the flow fields. Array 
Steven again tends to smooth and weaken the meridional fields. The actual error 
tends to be larger in the meridional estimate than in the zonal estimate. 

Figure 4-22 compares the actual zonal flux fields derived from array 


Harry. These estimates compare well with the actual values. Within the array, 
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point estimates of flux are good, but degrade slightly near the boundaries. 
Meridional flux estimates are also handled well by array Harry as can be seen in 
Figure 4-23. Array Harry provides an adequate resolution of the current fields, 
and thus the estimate of flux is improved . 

In summary we have compared the actual and estimated zonal and 
meridional components of heat flux at three model levels above 500 m. Point 
estimates of heat flux are good for the zonal case but very poor for the meridional 
case for array Steven. This is due primarily to the poor resolution of the meridional 
current field. We show a significant improvement in both zonal and meridional 


current estimates with array Harry due to its superior resolution. 
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Figure 4-19 Zonal Heat Flux: Comparison of the actual January 89 zonal 
heat flux and its estimate generated by array Steven. Units are in 


(watt meter-2) x 10 o 
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Figure 4-20 Zonal Heat Flux: Comparison of the actual January 89 zonal 
heat flux and its estimate generated by array Harry. Units are in 
(watt meter-2) x 10°. 
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Figure 4-21 Meridional Heat Flux: Comparison of the actual January 89 
meridional heat flux and its estimate generated by array Steven. 
Units are in (watt meter-2) x 10 is 
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Figure 4-22 Meridional Heat Flux: Comparison of the actual January 89 
meridional heat flux and its estimate generated by array Harry 


Units are in (watt meter-2) x 10°. 


ca 


2. Depth Integrated estimates 

The actual depth integrated heat flux for January 89 is shown in Figure 
4-23 and applies to both arrays. Large heat flux values are present in the western 
half of the array resulting from the warm, strong flow. Flow fields have been 
superimposed to show the direction of movement. The estimated depth- 
integrated flux field for array Steven has large errors especially near the center 
of the array as shown in Figure 4-24. The large positive error represents an 
under-estimate while the adjacent negative error represents an Over-estimate. 
The error is caused by the fine scale of the actual fields. 

In Figure 4-25 we present the January 89 depth-integrated heat flux 
estimated by array Harry. The estimate is superior owing to the array’s fine 
resolving power of the flow field. 

In summary we have presented a good and poor case of the 
tomographic estimate of heat flux in the ocean. The most significant error in the 
estimate of flux stems from the poor resolution of the current structure. Thus a 
tomographic system designed to measure heat flux must as a first priority provide 


good resolution of the flow field. 
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Figure 4-23 Depth Integrated Flux:Actual depth integrated January 89 heat 
flux for both arrays. 
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Figure 4-24 Depth Integrated Flux:Estimated depth integrated January 89 
heat flux estimated by array Steven. 
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Figure 4-25 Depth Integrated Flux: Estimated depth integrated January 89 
heat flux estimated by Harry. 
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3. Heat Flux Through a Transect 

We now present heat flux estimates through a transect in each array. To 
form these estimates, we take a zonal and meridional transect that bisects each 
array and then sum the depth integrated zonal and meridional components of heat 
flux along that transect. The transects of both arrays coincide so that we can 
make a direct comparison of the estimates. A number of transects are possible 
within the array, but we choose the center transect in each array to standardize 
our approach. 

Figure 4-26 displays the actual and estimated zonal heat flux through the 
center of the two arrays, Steven and Harry. In most locations along the transect 
we tend to under-estimate the zonal heat flux for array Steven. The estimates 
based on array Harry are far better. 

Figure 4-27 gives the actual and estimated meridional heat flux through 
the center of the two arrays. The estimates derived from array Steven are poor 
showing errors in both magnitude and phase. The estimates from array Harry are 
much better. 

Table 4.3 gives the net heat flux through transects in both arrays. Array 
Steven shows a 33% error in the zonal estimate while array Harry has a 4% 
error. Array Steven completely misses the meridional flux by being out of phase. 


Array Harry performs nicely with only a 9% error in the estimate of meridional 


flux. 


TABLE 4.3 Net Heat Flux: Comparison of the actual net heat flux and 


the estimates from both arrays in watts x 108. 


Component 
Zonal 
Actual 
Estimate 
Meridional 
Actual 


Estimate 


Steven 


32.0 
22.3 


5.0 
-4.8 
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Harry 


21.8 
20.8 


-10.2 
-9.3 


Zonal Heat Flux Watt/m’*x 10° 
Through a Transect 


Harry Estimate 
(J Steven Estimate 
BB Actual Fiux 


Zonal Flux watvm2 x 1068 


Figure 4-26 Net Zonal Flux: Net Zonal Heat Flux through a transect common 
to both arrays. 


84 


Meridional Heat Flux Watt/m?x 10° 
Through a Transect 


WE Actua! Fiux 
Steven Estimate 
Harry Estimate 


Meridional Heat Flux watt/m2 x 10e8 


Row Numbers 


Figure 4-27 Net Meridional Flux: Net Meridional Heat Flux through a 
transect common to both arrays. 
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Vv. CONCLUSIONS 


In this final Chapter we would like to offer a summary of our findings along 


with some conclusions. 


A. SUMMARY 
1. Covariance of the Solution 

One important finding is that although the covariances specified for the 
sound speed perturbation and flow fields are inexact the estimatés are only 
affected slightly providing that the system resolution is adequate. The actual error 
in the sound speed perturbation field was much less than predicted by the 
estimator. However, the current field estimates contain slightly larger errors than 
predicted by the estimator. This is due to the poor assumption we have made 
about the current fields themselves. Namely, the actual velocity fields are not 
isotropic. The estimator tends to redistribute the energy. As a consequence, the 
field is smoothed. Another problem is that we have assumed that the flow 
averaged 10 cm/sec at each depth when in fact the average current decreased 
with depth. In this case the estimator performed accordingly by under-estimating 
the near surface levels while over-estimating the deeper levels. 

Correcting these problems is a simple matter of specifying more 
accurate information to the estimator through the covariance of the parameter. 
For example we can specify a decrease of energy with depth for the current in 


addition to anisotrophy. 
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2 Acoustic Channel and Vertical Resolution 


The acoustic channel also plays a significant role in the application of 
tomography to a given region. In this simulation, we saw the vertical resolution 
degrade significantly for the August 89 profile. This affects the near surface 
estimate of heat flux. In summer months we might rely on a technique that 
couples tomography with satellite derived SST data to adequately resolve the 
surface temperature structure. 

3. Array Geometry and Horizontal Resolution 

Another important consideration is the selection of an appropriate array 
geometry. Horizontal resolution and array geometry are inter-dependent. In this 
study we found that for circular arrays, rather sparse transceiver spacing was 
adequate in resolving temperature. However, array size had to be reduced 
significantly to adequately resolve the current fields. This means that to deploy an 
array in a flow field with a complex structure we must consider reducing the area 
covered or increasing the number of transceivers. The the minimum spacing of 
transceivers is dictated by the spacing necessary to resolve current. 

An asymmetrical array was also simulated to determine what 
improvements if any could be realized from this type of an array. The 
asymmetrical array with the major axis oriented with the flow did provide a slighi 
increase in resolution. This increase in resolution, was small compared to the 
resolution increase from a reduction in array diameter. 

4. Comparison with Conventional Techniques 

Tomography compares favorably with conventional heat flux 

measurement techniques. We have demonstrated that the tomographic array is 


capable of adequately resolving both temperature and current fields over a large 
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area. We have shown that a six element tomographic array can adequately 


resolve heat flux in approximately 50,000 km? area of ocean. To provide similar 
coverage with direct techniques would require the deployment of 20 current 
meters to attain a similar resolution. 

Tomography also features the advantage of providing a noise filtered 
measurement of current. Conventional techniques are susceptible to 
contamination of data from secondary sources such as tides and internal wave 
motion. 

One disadvantage of tomography is that it is susceptible to degradation 
caused by failure of transceivers. Another words, failure of an array element 
could wipe out a large portion of data. We did not consider transceiver failure in 
this assessment although the topic is discussed by Kao (1989). 

5. Viability of the technique 

Tomography does offer a unique way of measuring the ocean’s 
temperature and flow structures. We have seen some of it’s strong points and 
weak points. Based on these preliminary results we have confidence in the 


technique and have specified ways of improving the estimate. 


B. CONCLUSIONS 
Acoustic tomography is a viable method of estimating ocean heat flux. It 
offers advantages of good spatial and temporal resolution. The method adequately 
resolved heat flux in this numerical simulation in the Norwegian Sea. The 
technique is flexible and should be adaptable to most ocean environments. 
Resolution of the temperature and current fields is the key issue in this 
experiment. A six element circular array, 250 km in diameter, adequately 


resolved both temperature and current over an area 50,000 km?. Heat flux 
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through a transect at the center of the array was estimated with a 4% error for 
the zonal component and a 10% error for the meridional component. 

A resolution analysis of the two circular arrays revealed that the minimum 
resolution length is roughly 20% of the array diameter. We also found that the two 
unknowns, &c and V have different variability in the ocean considered and require 
different minimum resolution lengths. The current field is the most difficult 
unknown to resolve and thus it is the limiting factor in the design of a tomographic 
array. 

Temperature fields were adequately resolved, even for a resolution length of 
120 km. Current could only be resolved adequately if the resolution length was 
less than 50 km. Errors in the current estimates related directly to system 
resolution. Point flux estimates were good as long as the current was well 
resolved. Depth integrated estimates show that vertical resolution was adequate 
_for the January 89 fields. The August 89 sound speed profile featured a poorer 
vertical resolution due to the channeling nature of the acoustic environment. 

The Gauss-Markoff estimator is only optimal when supplied with correct a 
priori information. Estimates of current were degraded by the poor assumption 
made about the nature of the current fields. Estimates can be improved by 
specifying an anisotropic current field having an energy level that decreases with 
depth. 

In conclusion, tomographic estimate of ocean heat flux is a technique that 


works and warrants further application. 
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